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ABSTRACT 



The recent failure of the scenario which uses only 10, 126 Higgs multiplets to 
yield large enough neutrino masses in the SO(10) MSGUT motivates an alternative 
scenario PQ where the 120 -plet collaborates with the 10 -plet to fit the dominant 
charged fermion masses. The small 126 -plet couplings give appreciable contribu- 
tions only to light charged fermion masses and enhance the Type I seesaw masses 
to viable values. We analyze the 2-3 generation core of the complete hierarchi- 
cal fermion mass system in the CP conserving approximation. Ansatz consistency 
requires — m s = m r — m M at the GUT scale Mx and predicts near maximal 
(PMNS) mixing in the leptonic sector for central values of charged fermion parameters 
and for wide ranges of the other relevant parameters : righthanded neutrino masses 
and relative strength of contributions of the two doublet pairs from the 120-plet to 
the effective MSSM Higgs pair. These features are preserved in the CP preserving 
3 generation system whose results are previewed : an additional consistency require- 
ment = #i2#23 011 the CKM angles at Mx arises for 3 generations. Right handed 
neutrino masses in such scenarios are all less than about 10 12 GeV. 



1 Introduction 

Renormalizable supersymmetric SO (10) GUTs[2, 3] have recentlyjl] been a focus of 
intense interest both as regards the GUT scale symmetry breaking 4 , spectrajSJ 
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[7| [H|,RG evolution[S] and matter Higgs Clebsches [SJ El- Even more interest arose 
due to the demonstration of the feasibility[9j of the generic Babu-Mohapatra(BM) 
program [Til] for a completely realistic fit of all the fermion mass and mixing data using 
only the 10 and 126-plet Fermion Mass (FM) Higgs. This scenario was initially also 
considered to be "predictive" in the neutrino sector. Detailed analysis [TT1 IT2*1 ITS] has 
shown, however, that successful charged fermion fits in fact require (arbitrary) choices 
of the (many) free phases present in the problem. Thus the successful fits can at most 
be considered indicative of likely MSSM C GUT embedding angles/ phases [Ti] when 
the leptonic mixing and measured neutrino mass splittings are taken as data. These 
embedding angles are physical in the sense that they affect [Hj the baryon violation 
and Lepton flavour violation rates in the Susy GUT and may therefore eventually be 
measurable. The renewed interest and excitement in the Babu-Mohapatra scenario 
began with the demonstration [9] that a natural linkage existed -at least regarding 
magnitude if not exact maximality - between Type II Seesaw JH] dominated neutrino 
masses and the manifest approximate b — r unification at scales ~ Mqvt observed 
in the MSSM and the near maximal neutrino (PMNS) mixing^Tj in the 2-3 sector. 
The nearly successful generalizations of the generic fits to the 3 generation (real) and 
realistic 3 generation (complex) case JJJ and the "completely" satisfactory generic fits 
when a small perturbation of the charged fermion masses by the 120-plet representa- 
tion was allowed [T2*] seemed to elevate this generic fitting scenario to near canonical 
status : particularly due to the natural ease of the Bajc Senjanovic Vissani b — r- Type 
II seesaw- large mixing linkage; even though the near exact maximality of the 2-3 sec- 
tor leptonic mixing does not have any obvious cause within such scenarios. The early 
observation pH HTi] that the Type I and Type II neutrino masses tend to be too small 
(and therefore |19j) may be considered to motivate more complicated GUT Higgs sec- 
tors) did not dampen the general enthusiasm. However, nearly simultaneously with 
the demonstration J2| that the Type II weak 10 © 126 ©weak 120-plet FM fit was 
completely satisfactory at least as far as neutrino mixing angles and mass squared 
splitting ratios are concerned, came the surprising finding [T3j that Type I and Type 
I © Type II fits are -in principle -equally viable. Furthermore we showed, nearly 
simultaneously [14], that Type II seesaw fits were most likely always far sub-dominant 
to Type I seesaw and Type I itself was too weak in the BM 10 scenario, at least in 
the fully specified and calculable context of the MSGUTj4J|H], and under an assump- 
tion of the genericity of the magnitude of a seesaw matrix derived from a successful 
Type I seesaw fit available in the literature [1311], even at the special points where the 
seesaws are weighted strongly [H EJ. Furthermore, with improved parametrization of 
the MSGUT symmetry breaking [20] we showed [JJ, by a detailed survey of the param- 
eter space of the MSGUT -assuming the typical values of the 10 — 126 parameters 
available in the literature [T2*[ ITS] were generic - that, not only was Type II seesaw 
subdominant to Type I Seesaw everywhere in the parameter space, but also Type I 
seesaw itself could not yield neutrino masses in excess of about .005eU. This failure 
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was traced PP to the double load borne by the 126 FM Higgs irrep and motivated us 
to propose the scenario that the 10-plet and 120-plet shoulder the main burden of 
the charged fermion fit, while the 126-plet coupling to the matter 16-plet is small 
so as to enhance its contribution to the Type I Seesaw (and lower that to Type II) 
without any strong (i.e involving 2-3 generation masses) lower bound on the mag- 
nitude of the 126-plet coupling due to the requirements of the charged fermion fit. 
Although the 120-plet has been considered-somewhat cursorily- previously in the 
literature [2JJ and has also recently attracted attention not only as a perturbation 
to correct "small defects" in the generic Babu-Mohapatra(BM) fits|12j but also as a 
dominant contributor to the charged fermion fits[22j in non-supersymmetric theories 
with radiative neutrino masses arising[23j from exchange of 16n-plet Higgs . Yet the 
combination of 10 — 120 — 126-plets with the particular- and experimentally well 
motivated - assignment of roles, in the context of the renormalizable Supersymmet- 
ric SO(10) GUT, that we proposed^ has not, to our knowledge, been previously 
considered in the literature. 

This is the first paper of a series in which we examine we examine this scenario ab 
initio and show that it enjoys a number of the same virtues manifest in the BM-Type 
II case [0] and also shown to exist in the BM-Type I caseU3j while being free of their 
overall neutrino mass debility. Moreover it leads to a restrictive type ofb — r = s — fi 
unification that radically reduces the latitude in choosing m^Mx) to the same level 
as the uncertainty in m T (Mx) yet proves not only compatible with la limits but also 
leads to an extremely robust prediction of very near to maximal (sin 2 9^ > .95) 2-3 
generation lepton mixing completely compatible with experiment and unlikely to be 
modified when first generation masses and CP violation are also accounted for! 

We first briefly summarize the essence of the reasons for the failure of the 10 — 126 
scenario |14[ IT] and give the fermion mass formulae in the case where a 120-plet is also 
present. We then analyze these formulae in the case where 10 — 120 plet completely 
dominate the charged fermion sector for the real two ( 2nd and 3d) generation case: 
where a completely explicit and analytic treatment is possible. This provides a insight 
into the the dominant core of the fermion Hierarchy and thus a clear paradigm for 
the analysis of realistic 3 generation case whose complications can be tackled by a 
perturbative treatment for which the 2 — 3 sector we analyze is a very stable, robust 
and non singular core and support which dominates the perturbative equations up to 
cubic order in the Fermion hierarchy parameter e ~ Q\ 2 ~ V^23 ~ 0.2. We formulate 
the fitting equations in a way adapted to determining their solution by an expansion 
in this e parameter. For dominant 2-3 sector neglect of the 126 couplings in the 
charged fermion sector is shown to yield an acceptable scenario with b — r = s — fi 
unification and maximal neutrino mixing and masses for precisely the expected values 
of the 2-3 generation masses at Mx and a wide and plausible range of the remaining 
parameters (M^ and GUT doublet/MSM doublet Higgs fractions) ! We conclude 
with an preview of the results of paper II where the 3 generation real case is analyzed 
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using the perturbation method based on our results here. 

2 Difficulties of Babu-Mohapatra Seesaw in the 
MSGUT 

In (TJ we wrote the Type I and Type II seesaw^H] Majorana masses of the light 
neutrinos as : 

Ml = (1.70 x 10~ 3 eV) nsinp 

Ml 1 = (1.70 x 10- 3 eV) F n /sin/? 

n = (h-mr'Ch-Sf) 

where v = YlAGeV , h, f are the Yukawa coupling matrices of 10, 126 to the 16 
plets containing fermion families and (3 is the MSSM Higgs doublet mixing angle. 
The functions Fj, F n are defined (upto irrelevant phases ) as 



10" Ax 19 , , [Z2 (1 ~ 3s) g' 3 
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\P2P3P5\ 



(x 2 + l) (4s-l)g 2 
rfX {x - 1) V z 1G \\ x(l - 3x) q' 3 q2P5 



Here 7, A, T], g are couplings ~ 1 while the Pi,qi, z% 
are certain polynomials |H 1^0*1 IT] in the variable x : which is the only "fast" parameter 
controlling the GUT scale symmetry breaking in the MSGUT[4, 6 . This MSGUT 
is based on the 210 — 10 — 126 — 126 Higgs system and is tuned to keep one pair 
of Electroweak doublets light in the MSSM that emerges as the low energy effective 
theory. The modification( Ax = Log\o(Mx / GeV) — 16.25 ) of the one loop unification 
scale due to GUT scale threshold corrections |Sj is constrained by the proton lifetime 
to be greater than — 1. 

The 'proof of |T3J proceeds by first observing that in typical BM-Type II 
fits|ll[ IT2"| ITH] the maximal value of / eigenvalues is ~ 10 -2 while the corresponding 
values for h are about 10 2 times larger. As a result h ~ 10 2 /. This implies that 
R = Fj/Fjj < 10~ 3 for the pure BM-Type II not to be overwhelmed by the BM-Type 
I values it implies. Furthermore in the BM-Type I fit of ^1 the requirement of large 
neutrino mixing yields, typically, h max ~ 5f max ~ .3 for the maximal eigenvalues of 
h . Thus hmax is much smaller than the naive estimate (^33 — f-^) 2 / f^z ~ 10 1 so that 
magnitudes Fj > 10 2 are required for neutrino masses ~ .05 eV. This deserves further 
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examination on the basis of either a catalog of Type I fits of the type found in[T5] 
or an analytic perturbative treatment of the fitting problem |24[ I33j. However, the 
reason for the lowering of the eigenvalues of n below the naive values appears to be 
the necessity to have roughly equal elements in the 2-3 sector to allow large mixing 
and so should be robust. 

The complete control we have over the MSGUTjU E21 HI] allows us to 
demonstrate |14[ P that both these conditions are not achievable anywhere over the 
complex x plane except where they also violate some aspect of successful unification 
(typically the requirement that Ax > —I when combined with the above requirements 
of the completely determined coefficient functions of the MSGUT ensures the failure 
to generate adequately large neutrino masses : see [I] for details). 

3 2 Generation FM fits in NMSGUTs 

The above recounted failure of the Babu-Mohapatra scenario, at least in generic 
cases of the MSGUT, could provoke one[in| to entertain more complicated scenarios 
which give more freedom to choose the parameters that enter the neutrino mass 
formulae. The likelihood of Type II failure, already noticed in ^H] there motivated 
the introduction of additional 54 Higgs as a work-around of the obstruction posed. 
However, as we already noted in a proper demonstration would then require a 
complete recalculation, of the various mass matrices and coefficient functions which is 
-and may remain- unavailable, not to speak of the new essentially 'fast" parameters 
introduced by the modified -and considerably more complex GUT scale SSB scenario 
of such work-arounds. To us, allowing in the arbitrarily excluded 120-plet, particularly 
since it does not destroy the hard-won solution of the GUT scale SSB problem, is far 
more palatable and cogent. 

If this is done the main effect is to introduce two additional doublet pairs, from 
the (1, 2, 2)and (15, 2, 2) Pati-Salam submltiplets of the 120-plet. The GUT scale SSB 
is undisturbed. The Dirac masses in such GUTs are then generically given byjHJG] 

m u = v(h + f + g) 

m u = v(h-3f + (r 5 -3)g) =v(h-3f + r' 5 g) 

m d = v(rxh + r 2 f + r 6 g) (2) 

m l = v{rih — 3r 2 f + (f 5 — 3f 6 )g) 

The form0E] of the new generic coefficients r^,^^ in the particular case of the 
NMSGUT is of some interest for future work|18 : 



rs = — 7 • Fx ; r5 = — 7 • m cot/3 

o.Q + iyia^ «6 + iv oat*, 
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r 6 = — , 7= 5 cot/3; <? = 2icn/ -(a 6 + zv 7 ^) sin/3 (3) 

a 6 + zVoa 5 V 

here g^s is the 120 — 16 — 16 coupling,and the are constants determined in 
terms of the the cei(azi) which are are the fractions of the Electroweak doublet H[l, 2, 1] 
contributed by the two [1,2,1] ( [1,2,-1]) doublets in the 120-plet . We have omitted 
similar details for the 'old' coefficients r», f»;i = 1,2[S]. Note that with free our 
analysis is generic and applies to any 50(10) theory with the 10 — 120 — 126 FM 
Higgs system. 

The right handed neutrino mass is M p = fa and the type I seesaw formula is 

Ml = vr,h ■ h = (h-3f)f- X (h-3f) (4) 
where a = ta ^ a and a is the GUT scale vev of the 126. 

c*2 sin p 

The essence of our proposal is to neglect the presence of / in the Dirac masses of 
the 2-3 generations on the basis of an assumption that h,g » f. Notice that this 
implies the strengthening of Type I at the expense of Type II Seesaw mechanism. 
When one analyzes the complete 3 generation case one finds that dominance of h, g 
can only be strong enough to justify the neglect of / to 0(e 3 ). At 0(e 4 ) non zero 
values of / are required in order to preserve the consistency of the charged fermion 
mass fit at and beyond fourth order. Note that such small / values imply that the 
typical right handed neutrino masses would lie in the range below 10 12 GeV. Since 
these masses are also thought to obey the Davidson-Ibarra bound M VR based on 
reheating constraints in the popular leptogenesis scenario [26 for baryon asymmetry 
generation (which is so natural in theories supporting Seesaw neutrino masses). This 
may be a useful clue provided by renormalizable MSGUTs towards constraining the 
parameter space of the leptogenesis scenarios based on the Type I seesaw mechanism. 
In any case in the present 2-3 generation Real Core analysis we will set / to zero and 
take all parameters and Unitary matrices to be real. The required generalizations 
to the 3 generation case are quite easy -although computationally tedious - and are 
given in the next paper of this series j23J. 

The mass terms above must be matched to the renormalized mass matrices of 
the MSSM evaluated at the GUT scale. While doing so one must allow^I] for the 
possibility that the fields of the GUT are only unitarily related to those of the MSSM 
at Mx- This introduces several 3x3 unitary matrices as well as various phase 
matrices into the fitting problem. Thus when matching the charged fermion dirac 
mass matrices (in which we neglect the contribution of thel26) we get 



m u = v{h + g) = V^D u Q 

m d = v{nh + r G g)VjD d R 

m 1 = v(rih + r 7 g) = D\L 

r 7 = f 5 - 3f 6 (5) 
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Where D U)dd are the charged fermion masses at Mx and V u , Q,Vd, R = C^Q, L, Vi 
are arbitrary unitary matrices (C is the CKM matrix) which are to be fixed by 
convention -where allowed by conventional ambiguities - or determined by the fitting 
procedure in terms of the low energy data, or left as parameters to be determined by 
future experiments sensitive to degrees of freedom and couplings (e.g baryon violating 
couplings) that the low energy data is not. 

To put our equations in a form transparent enough to clearly separate the con- 
tributions of the 10, 120 (and -eventually- 126) plets,24j we write the matrices V Ut d,i 
associated with the anti-fermion fields as new unitary matrices & u ,d,i times the fermion 
field matrices Q,R,L. 

V d = <$> d R ■ V U = $ U Q ■ V t = ^L (6) 

We then separate symmetric and antisymmetric parts for each charged fermion equa- 
tion 

Z = <$> T D + D<& ; A = $ T D- D$ (7) 

Then it is easy to solve for the matrices h,g using (say) just the d-type quark 
equations and substitute in the other two equations to obtain the equations in the 
form 

C*Z d & = n Z u ; Z d = DZ l D T ; D = R*L T (8) 
C*A d tf = r 6 A u ; r 7 A d = r,DA l D T (9) 

In the 2 generation case the antisymmetric equations© serve only to fix the 
parameters r^^r-j and thus play no further role. On the other hand since C,D are 
orthogonal matrices it follows that r% = Tr(Z d ) /Tr(Z u ) so that we can write these 
two equations in the (dimensionless) form 



Si 



CZ d C T Z u 
TrZ d TrZ u 



- = Z d -DZ,D- = 

TrZ d 1 ' 

The convenience of this form is that the hierarchy within each generation will now 
structure the problem and render its solution almost trivial in the two generation 
case and tractable in the 3 generation case. It is amenable to an expansion in the 
hierarchy parameter e which reduces the problem of finding a sufficiently accurate fit 
to the ~ 10 4 intra-generational variation in mass to an exercise in going to sufficiently 
high order in the expansion parameter e[~ 
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3.1 Analytic Solution of the Charged Fermion Mass Rela- 
tions 



We parametrize the matrices $ n ,d,i, C,T> as 



cos Xu,d,i sm x u ,d,i 
- sin Xu,d,i cos Xu,d,i 



CV = I COSXc ' v sm Xc,v \ (11) 
\ -siriXc.D cosXc,v J 

and define Am = u 3 — u 2 , T u = u 3 + u 2 , Ad = d 3 — d 2 , T d = d 3 + d 2 , Al = 
h — h , Ti = l 3 + l 2 . 

Then the 2 independent equations in S\ = reduce to 

t&riXu = tan(xd - 2x c ) 

cosxd = acos(xd-2xc) (12) 
(Ad) M (d 3 -d 2 ){u 3 + u 2 ) 



a 



(T d )(Au) (d 3 + d 2 )(d 3 -u 2 ) 



Where Xc is the Cabbibo angle in the 2-3 sector i.e in the toy CKM matrix C and 
^2,3, d 2} 3, l 2 ,3 are the charged fermion masses of the 2nd and 3d generations up to 
arbitrary signs. Since Xu = Xd~ 2>Xc and Xc ~ 0(e 2 ), Xud ~ Xd to order in e 2 . The 
second equation in (|T3j) equation can be solved for tan Xd '■ 

esc 2xc 

tanxd = cot2x c (13) 

a 

Since a = 1 + (e 2 ), it is clear that the leading (~ e~ 2 ) contributions cancel leaving 
behind an 0(1) result. Similarly the three equations S 2 = yield 

tan(x; - 2xd) = tan Xd 

Td d 3 + d 2 

cos xi = cos Xd = , , , cos Xd 
J- 1 h + h 



(14) 



Thus 



Xd = (Xi-Xd)/ 2 ; Xz = ±X« 

_ 1 T d cosXd M _v 

X« = cos (15) 

J-i 

To leading order xi — However note at at next to leading order the ratio Tj/Ti 
must not become so large (tan%° ~ .9) as to make Xi the inverse cosine of a number 
bigger than 1. 
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The third equation yields the all important consistency condition 



d 3 - d 2 = Ad = Al = l 3 - l 2 



(16) 



i.e 



m b (M x ) - m T (M x ) 



m s (M x ) - m M (M x ) 



(17) 



This is one of the main results of this paper. It should be understood that there is a 
sign/phase ambiguity-arising from the arbitrariness in the unknown Unitary matrices 
that relate the SGUT and MSSM fields (see above)- for each of the masses, while 
their magnitudes are determined by the RG flow of the MSSM into the UV. 

4 b — r = s — /i Unification 

Several questions arise regarding the status of this version of b — r unification which is 
imposed by the consistency of the assumption that the 10, 120 multiplets dominate 
the heavy fermion mass matrices. Firstly it should be noted that the necessity to 
the zeroth order(in rj = e 2 ) version of this relation i.e m^Mx) — m T (Mx) = was 
already noticed by Bajc and Senjanovic 22j while studying a technically related model 
but with a but quite different [non- super -symmetric) scenario. In that model the 
10 — 120 are used to fit the 2-3 generation masses and neutrino masses are generated 
radiatively [23] due to exchanges involving a a SO(10) spinorial Higgs 16#-plet. Thus 
the motivation, origin and properties of that model are totally distinct from our 
proposal which is supersymmetric (and thus protected from radiative corrections of 
its Superpotential) , employs no Higgs 16-plet and has a tree level Type I seesaw using 
the 126 (which also necessarily contributes to the fermion masses at order e 4 [21]) in 
the thxee(but not two) generation case. We have here analyzed completely the two 
generation case and to order e 6 for the 3 generation (CP conserving case ) in [2*lj . 

The first question that arises is whether this relation will be preserved when the 
first generation and then CKM CP violation are also introduced. The explicit ana- 
lytic solution trivially derivable from the above equations given above coincides, order 
by order in perturbation theory , including the b — r = s — ii constraint , with the 
solution found [21] by expanding all angles \ in powers of e and solving the equations 
Si = = S*2 order by order in e (and truncating all expansions by dropping undeter- 
mined coefficients left after solving to desired order). The b — r = s — ii unification 
constraint arises at order 77 = e 2 . In the three generation case the analytic solution 
is not available but the complete regularity observed in the perturbation theory in e 
and the extreme smallness of the first generation perturbations makes an expansion 
in e well motivated. Then we find that the b — r = s — 11 unification constraint 
is reproduced at 0(e 2 ). One may also expect that this constraint will persist when 
the CP violation is introduced and a perturbation theory in e used to solve the the 
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fermion mass relations since it is well known that CP violating effects are Wolfenstein 
suppressed. In that case, however, the relation will also contain phase uncertainties 
generalizing the sign uncertainties present in the real case. Thus b — r = s — fi uni- 
fication may be considered as one of the signature "tunes" of the fermion mass (FM) 
fit in the NMSGUT scenario we proposed^] and have begun elaborating on here. 

Moreover the pursuit of the perturbation theory to 0(e 3 ) gives [2"f] an additional 
remarkable constraint : 

#13 = #12 #23 ( 18 ) 

between the CKM rotation angles ! No further parameter constraint was found up 
to 5th order in e. However at 4th order we found that the contribution of the 126- 
plet becomes necessary to avoid inconsistency in the Si,S 2 fitting equations, which 
remain combined in a single equation aSi + bS 2 = rather than vanishing separately, 
while the coupling of the 126-plet is determined to be ~ S 2 ~ e 4 . Given enough 
computer power we see no reason why the perturbation cannot be carried to sixth or 
even seventh order in e : which should be more than sufficient to ensure convergence 
to the accuracy to which the lepton masses at Mx are known : given our uncertainty 
regarding the crucial quantities Ms, tan/3. 

Since d 2 /l 3 ,l 2 /h,Xc ~ and u 2 /u 3 ~ e 4 it immediately follows that the solution 
of eqn (jl3H13]) to leading order in epsilon is 

tan X °=tan X ° = ^ + 0(e 2 ) (19) 

Thus the magnitude of tan \d is determined to lie (see data values below) in the range 

|tan X °| = .4974(+.19)(-.l) (20) 

Thus the accuracy of the leading order is ~ 10%, while the experimental uncertainties 
are larger. For all practical purposes we may eliminate £3 in favour of tan%° as the 
basic parameter and let it vary in the range above with either positive or negative 
sign. Note that these values will prove crucial in showing that the leptonic mixing is 
near maximal over a wide range of parameters. 

4.1 Numerical Values 

Even the very precisely known low energy lepton masses become smudged by uncer- 
tainties when RG evolved past the unknown thresholds of order Ms and M p with 
unknown tan/?. Since the effects of CP violation on the 0(e 3 ) relation 6^ 3 = 61 2 @23 
may well be non-negligible, the precise numerical implementation of the constraints 
seems premature. Still a survey of the typical fermion mass-mixing data sets is quite 
comforting since it exhibits at least approximate compatibility with the constraints 
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found by us and makes clear the huge parameter space available for satisfying them. 
Moreover implementation of the necessary b — r = s — \i constraint implies a radical 
(10 — 100 fold !) reduction in the uncertainty allowed in rrib(Mx) due to its strict 
correlation with m T (Mx) '■ which is much more precisely known(see below). 

In j2Zj the following values are given for the two loop renormalized running fermion 
masses (in units of GeV) at /i = M x = 2 x 10 16 GeV for the representative high tan/3 
value tan (3(M S = ITeV) = 55 

m u = (0.72441S) x lO" 3 ; m c = .2105049i;K£ ; m t = 95.1486^;lf 

m d = (1.4967+ 2 ' 2 4 7 1 8 57 ) x 10~ 3 ; m s = (29.8135^96?) x 10~ 3 ; m b = 1.4167^™ 

m e = (0.3565;°;X) x 10~ 3 ; m„ = (75.2938;°;^) x i -3 . mr = 1.6292+°;^ 3 

tan/3 = 52.0738;™; ^ = 117.7947;™; Vd = 2.2620;™ 



It is apparent that after including the sign ambiguity |d 2 — h\ can be anywhere 
from .046GeV to .lOSGeV" with an notional accuracy limited by the low energy error 
induced uncertainty in m s (Mx) of about .OOAGeV. On the other hand \d$ — h\ is 
about .2 GeV±.5 GeV (!) (evidently ds,l 3 must have the same sign !). There is a 
long way to go before the validity of the constraint can be verified, even granted that 
no additional ambiguities due to CP violating phases arise. 

Our approach is thus to take as reference central values the relatively accurate 
values of m^ TtS , at large tan[3 (obviously favoured by SO(10) GUTs) (say tan/3 = 55 
for definiteness) and assume d 3 = m r ± m s ± m M w 1.63 + sign[d 2 ].03 — sign[l 2 ].075) 
to satisfy the consistency constraint, secure in in the knowledge that the data allows 
nib to lie in the range (1.23, 1.9) for tan/3 = 55 and somewhat lower values for lower 
tan j3. However we find that sign[d2\ = +, Sign[l 2 ] = — is impermissible since it leads 
to an imaginary value for xi when one solves eqn ()14|) . 

The values of d% required in our scenario by the combination of 6 — r = s — fi 
unification and near maximal leptonic mixing tend (at d% ~ 1.585) to be somewhat 
larger than the central value quoted for nib above but still comfortably within the 
la error bars. This can be considered a prediction of this model which may become 
constraining as limits improve. Note particularly that the implication of this type of 
novel constrained unification leaves only the relatively tiny la uncertainty (~ .001 — 
.OSGeV ) in m T (M x ) as a fudge factor while d 3 whose (i.e m(,(M x )'s) uncertainty 
could have been 10-100 times larger is tied down to vary in tandem with m T (Mx)- 
Yet we will see that the prediction of the maximal mixing is still extremely robust. 

The values of the mixing angles at Mx given by (who assume 5c km = 7t/2 
for convenience ) are Q\ 2 = .221,6* 2 3 = .037, #i 3 = .003, with no uncertainty quoted. 
The product 9i 2 9 23 = .008 ± .002 is of the same order of magnitude as 6*13. Since 
this relation arises in the 3 generation case at 0(e 3 ), by which stage the CP violation 
also enters it seems legitimate to hope for a considerable (~ 100%) role for Sckm in 
modifying the relation 6^3 = #12^23 between angles. 
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5 Neutrino Masses and Mixing 



In this section we will first determine the angle in the PMNS JHj mixing matrix and 
then proceed to show that this mixing angle can can easily satisfy the near maximality 
constraints ^7] imposed by experiment. 



5.1 PMNS Matrix 

Under the assumption that the couplings h, g of the 10 and 120 completely dominate 
coupling / of the 126, the Type I seesaw masses in this theory are given by 




f 

= r i R T YjRS- l Diag{l 1 p)S~ 1T R T Y d R 
= R T FR 

= L T VD U V T L (21) 

Where V is the Lepton mixing (PMNS) matrix in the basis with diagonal leptonic 
charged current, D u the light neutrino masses extrapolated to Mx, and Y d a con- 
venient dimensionless form of the linear combination of Z^, A d that determines the 
PMNS mixing. 

The explicit form of Y d is 

/ %d cos Xd - sin Xd - r' 5 j^ sin Xu \ 
Y d = (22) 

V -sinxd + r^sinxu 2(l + ^)cosxd/ 

So far we have not used the family basis ambiguity of the SO(10) GUT that allows 
us to perform an arbitrary unitary redefinition of the matter 16-plets at will. Since / 
is symmetric it may be written as / = S T D jS where S is unitary and Df is diagonal 
and real. The basis ambiguity can be fixed by a choice of S to be any given unitary 
matrix : for example to be unity. Here we make the convenient choice S = R since this 
removes the obscuring factors RS' 1 in the Type I mass formula eqn()21|) above, leaving 
behind the two eigenvalues of / as free parameters. One of these can be extracted 
into the overall scale say leaving behind the parameter p = f 2 / f% = M v /M p to 
affect the determination of the mixing angle xv- 

Imposing the family basis choice and diagonalizing the matrix F we immediately 
obtain the neutrino mixing and masses to be 

F = Y d T Diag(l,p)Y d = ^Dtag(F 2 ,F 3 )^ T 

V = D*T ; D u = ^Diag(F 2 ,F 3 ) (23) 
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The mixing angles x$ °f the 2-d rotation matrix T which diagonalizes F to 
Diag(F T , F±) are given by 



Xf = tan 1 (uj f ^Jl + u 2 



F) 



The eigenvalues F ± obey 



WF = fil^E (24) 

^12 



1 



F± = -{Tr{F)±^tr{FY-Det[F\) 



sign[Tr{F)] = signal - F 2 )] (25) 
Note that since (to leading order in e) 

Tr(F) = sin 2 X d(4 + 16p + r 2 (l + p) - 4r 5 (l + 2p) + 4pcot 2 xrf) (26) 
the hierarchy will invert for 

P (r 5 -4)2 + 4cot[ Xd ] 2 - P ^ 



5.2 Maximal 2-3 PMNS Mixing 

We now have a formula for the leptonic mixing in terms of two quite arbitrary pa- 
rameters p, r 5 = r' b + 3. The value p = 1 corresponds to completely degenerate 
right handed (heavy) neutrinos. While \p\ << 1 and \p\ >> 1 correspond to hier- 
archical righthanded neutrino masses, rs = corresponds to negligible contribution 
from the SU(4) singlet pair of MSSM doublets from the 120-plet (i.e Pati-Salam sub- 
representation representation (1, 2, 2) C 120) to the low energy symmetry breaking. 
Thus it may be realizing a Georgi-Jarlskog mechanism|2^j for the 10 © 120 based fit 
of the 2-3 generation charged fermion masses. We will show that 

• For any value of r 5 there exists a range of values of p for which the the leptonic 
mixing is near maximal (1 > sin 2 2\p > .9). 

• There exist ranges of r 5 where the width of the p-band where 1 > sin 2 2xp > -9 
becomes very large. 

Thus large mixing is always achievable and inevitable in certain broad regions of the 
r 5 , p plane. 

From the analytic solution of the charged fermion fit feqns (jl3H4j) and eqn (}2"2"|) we 
find that the two possible values of the total mixing are 



xF = x~ F ; x { P } = XF + Xd (28) 

13 



So that the two possible values of the mixing parameter sin 2 2%p (which is indifferent 
to inversion of the vl hierarchy ) are 



sin^xg = 



1 



sin 2 2x? = 



(1+^) 

(- 1 + 2u F tan xa + tan 2 Xd) 
(l+cc; 2 )(l + tan 2 Xd) 2 



(29) 



sin 2 2x^ is maximal at cup = uj- 



(0) 



0, while sin 2 2x^ is maximal at 
2tanxd 




2 tan 2 Xd - 1 



(30) 



Moreover the range of ujp values which leads to near maximal mixing (1 > sin 2 2x P > 
.9) is u F E (-1/3, 1/3). While for 1 > sin 2 2xp > -9 this range lies between 



There are in principle different combinations of signs : sgn{d 2 , l 2 , u 3 , u 2 , Xc} that can 
be chosen for the various parameters. They make no difference to the mixing at the 
leading order. However, the combination d 2 positive, l 2 negative, leads to imaginary 
Xi in the full theory(i,e beyond leading order where xi — ^Xd) and hence must be 
discarded. The remaining sign choices give results that are only marginally different 
from those where we take all masses and Xc as positive. So will not here discuss their 
minor differences and focus only on the "all positive" sign choice in the interests of 
clarity. 

5.3 Depiction of the Large Mixing Parameter Regions 

The regions of the r^,p plane that support large PMNS mixing can be fairly com- 
pletely and accurately delineated by working to leading order in e. Then tanxd ~ 
ms/xc^r is the only input parameter required. In leading (and determinative) order 
the formulae are not obscured by the other parameters whose role in determining the 
PMNS mixing is quite marginal. Note, however, that since we have an exact analytic 
solution numerical work can always use it to confirm and refine the features visible 
at leading order. 

It is easy to show that to leading order in e 



Note particularly the singularities at p = and at r 5 = 4 and the fact that the 
result is dependent only on r 5 ,p, tan Xd- There is clearly never a large mixing Xp"* 



(1,±,.9) 



±t 2 + 6tT 1 
3t 2 ^2t-3 



(31) 



(-(-2 + r 5 ) 2 + (-4 + r 5 ) 2 p - 4 p cot(xd) 2 ) tanfrd) 
4 (-4 + r 5 ) p 



(32) 
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solution for p = or r§ = 4, and it is easy to check that for these values the limit as 
uj — > oo, sin 2 2x^p — > which never yields large mixing for tanxd G (.35, .6) 

which is the experimentally allowed range. Clearly eqn ()33|) implies that for arbitrary 
(r 5 , tan Xd) one can always find a value of where the mixing sin 2 2xp* is maximal. 

Pma * (-4 + r5) 2 -4cot( X ,) 2 [ ' 

For example, when r 5 = 0, tan^d = .4974 , p^ ax = -23.8511 (while p inv = -0.124348, 
so the hierarchy is inverted but not acutely so). Moreover even for p = l,r 5 = 
(degenerate ur, Georgi-Jarlskog point) one finds 

uj f = - (cot Xd ~ 3 tan Xd) (34) 

which gives sin 2 2x^p = .9835 at tanx rf = d 2 / (hXc) = -4974. Giving already a hint of 
the robustness of large mixing. 

For given rs, t&nxd the edges p(°> ±1 /3) of the large mixing band for the mixing Xp^ 

are 

„(o,+,.9) _ 3 (-2 + r5) 2 tan(xd) 



-4 (-4 + r5 + 3 cot(xd)) + 3 (-4 + r5f tan(x. 



(0.-..9) 3(-2 + r5) 2 tan(xd) , . 

-12 cot(xd) + (-4 + r5) (4 + 3 (-4 + r5) tan(xd)) 

while the width of the band is 



A (0,.9) = p (0,|) _ p (0-|) 



24 (-4 + r5) (-2 + r5) 2 tan(xd) 



144 cot(xd) + (-4 + r5) 2 (-88 + 9 (-4 + r5) 2 tan(x 



(36) 

An exactly parallel discussion can be given for the other case, we omit the tedious 
details. 

The (simple) poles of A^ - 9 ) are at r5 -> -1.57871,1.10189,6.89811,9.57871. In 
Fig.l we plot the width of the large mixing stripe versus r 5 to illustrate its large 
variation. It is apparent that it is precisely the poles of A^ 0v9 ^ that cause the width to 
become large and that the whole region between r% = —1.6 and r§ = 10 (particularly 
r 5 G (6.5, 10) supports very large regions where the mixing is maximal. 

This is clearly seen in the corresponding contour plot Fig. 2 where we have 
truncated the extent in the p direction to avoid obscuring the small mixing band 
around p = 0. 

An exactly parallel discussion can be given for the mixing Sin 2 2xp^ , we omit the 
details and give only the contour plot in Fig. 3, which has an analogous structure to 
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Figure 1: Plot of A^ '- 9 ) vs r 5 , at lowest order in e, at fixed tanxd = .4974. The poles 
of A^ '- 9 ) cause regions of a very wide range of p to support large PMNS mixing. 



that for Sin 2 2xp ■ Note that the the very narrow width of the large mixing band for 
r 5 G (1.3,5) has caused it to become invisible due to the large range on the p axis. 
It is easy to see this relatively narrow band by magnifying the plot. In Figure 4 we 
show an example of the variation of the mixing parameter with tanxd keeping r 5 ,p 
fixed at values in the maximal mixing regions (sin 2 2x~p > -95). 

Although the exact formulae given above have a completely regular expansion in 
the small parameter 77 = .04 it is of interest to see how much the leading order results 
are modified when we include all orders in 77 and also to examine the effects of varying 
the other charged fermion masses e.g m top . The explicit formulae -although trivial 
to write down- are not very transparent . However plots (Figs 5-7) analogous to the 
ones for the leading order formulae tell almost the the whole story : 

• The contour plots for the exact results are essentially identical to those for 
the leading order case apart from minor changes in fine structure and a slight 
squeezing of the maximal mixing regions. The similarity is so great that we 
show only one example as Fig. 5 : the contour plot Sin 2 2xp\ We show the 
plot for a different sign choice for d 2 , h to emphasize that all these (sign choice 
and higher order effects) cause very minor modifications as is apparent when 
one compares Fig. 3 and Fig. 5. 
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Rho 



Tanchi_d = .4974 




r5 



-5 5 10 15 

Figure 2: Contour Plot of the mixing parameter Sm 2 2#23 = Sin 2 2xp\ to lowest 
order in e, on the (r 5 , p) plane at fixed tan^d- The contours shown are at Sin 2 2xp* = 
.85, .9, .95 with the shading becoming darker towards smaller values. 

• The variation with m t0 p{Mx) is quite marginal. An example of how it affects 
the mixing near a maximal mixing point is shown in Fig. 6. 

• Similarly the variation with \c near a maximal mixing point is shown in Fig. 7. 

• The role of m c is essentially negligible since it is 0(e 4 ) relative to m t and the 
equations are sensitive only to their ratio. 
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Figure 3: Contour Plot of the mixing parameter Sin 2 2023 = Sin 2 2xp' , to lowest 



order in e, on the (r 5 , p) plane at fixed tanx^. The contours shown are at Sin 2 2x^p 
.85, .9, .95 with the shading becoming darker towards smaller values. 



6 Discussion, Conclusions and Outlook 



The work presented in this paper is based upon the very productive insight [H] that 
it is the mixing in the 2-3 sector of the Fermion mass matrices that is the true stable, 
organizing core of of the entire complex hierarchical fermion mass matrix. Possibly 
due to the fact that the mixing angle #23 m the CKM matrix is much smaller than the 
Cabbibo angle 8^ 2 this realization has been slow to be digested by the many workers 
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{ rho^l , r 5^ } 




fascinated by the the problem of finding a non-trivial insight into the apparently finely 
crafted inner workings of this wondrous "primal artefact" . 

Counter to the fact that #23 << $12 one m ust consider the extreme tininess of the 
first generation masses relative to those of the 2-3 generations, which ensures that we 
must begin by understanding the smaller mixing angle. It was not until the discovery 
of the tiniest mass (differences) , via the observation of neutrino oscillations, precisely 
in the range expected from a seesaw between the Electro- Weak and GUT scales that 
the deep inner connections between the physics of the Largest (GUT), Smallest (m„), 
and "Everyday Electroweak " (Mw) mass scales emerged into plain viewjSU EH E] 
The initial surprise at the natural explanations lurking in the analysis of the simple 2- 
3 system and the apparently unique and natural simplicity P of the Type-II Seesaw, 
were followed by further surprises at the demonstration that similar large mixing 
could[13j also be obtained ("sporadically"?) in some of the possible solutions of the 
generic BM[10 j fitting problem in Type I seesaw (or a combination of both seesaws). 
The price paid for this "numerical discovery" was that it was based on tuning the 
many unknown phases present in these models to achieve the desired large neutrino 
mixing. 

Unfortunately the BM program appears to meet [HI EJ |2H 120] obstacles to its 
yield of large enough neutrino masses both in its Type I [3J E] and Type II 
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{u2->.21,u3->95, d2->- . 03, d3^1 .525, 12->- .075, 13->1 . 63 , chic^ . 037 } 




Figure 5: Contour plot of the exact solution for the mixing parameter Sm 2 2#23 = 
Sin 2 2xp on the r§, p plane, at central values of the charged fermion parameters (given 
on plot). The contours shown are at Sin 2 2xp = .85, .9, .95 with the shading becoming 
darker towards smaller values. The white region corresponds to sin 2 2\p > 0.95. Note 
the {sgnfa], sgnfa]} = { } sign choice. This may be compared with Fig. 3. 



[THl CH HI EHl ED] versions at least in the most desirable context of the MSGUT. 
Still this experience with the fitting problem provided the valuable lesson that the 
phase complications of the full complex CKM matrix problem would respect in large 
part the pattern discovered in the real two and three generation fitting exercise [T3] 
: which expectation is one of the operational assumptions of our work. 
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■ {u2->.211, 13-»1 . 63, d2^. 03, 12-» . 7 5 , chic^ . 37 , rho^2 .3, r5-»-l} 




m_top 



Figure 6: Plot of the exact solution for the mixing parameter Sm 2 2#23 = Sin 2 2xp , 
vs m top (Mx) , with sign choice {sgn[d^[, sgn]^-^) = {++} 

The likely smallness of the seesaw masses was noticed long before our work [3J E] 
and was considered in ^!|| to be a motivation for constructing more elaborate GUT 
scale Higgs structures to evade the tight constraints jH |H1 EH EE] imposed on seesaw 
coefficients by the very simple yet fully functional (and fully calculable jSJ E] !) 
structure of the MSGUT. Our contribution was a systematic survey of the entire 
parameter space of the MSGUT to show that the Type II contribution could not 
be saved by enhancement at the special points of the MSGUT parameter space and 
that the same also went for Type I. These conclusions were based on the genericity 
of the values extracted from Type I, II fits available in the literature and need to 
be confirmed on a wider sample generic fits or more satisfactorily by an analytic 
(perturbative) solution of the fitting problem in the MSGUT 24J. We are personally 
reluctant to enthuse about the assumption of more complex AM Higgs (i.e GUT scale 
SSB ) structures. A "shotgun " programmatic attitude towards decoding the enigmas 
of unification seems questionable, because if there is no tight discipline restraining 
"multiplication of hypotheses" (representations) while speculating on extensions of 
established theories then it is difficult, if not impossible, to know what one would be 
testing, if one ever could begin to do so : the wideness of the space of extensions or any 
specific unavoidable implication of theory. Since our knowledge about nature is sure to 
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Figure 7: Plot of the exact solution for the mixing parameter Sm 2 2#23 = Sin 2 2xp , 
to lowest order in e, vs chi c with sign choice {sgnfa}, sgn^]} = {++} 



arrive only incrementally and after hard struggle to define and test every additional 
observable, it should be kept in mind that extensions beyond those "minimally" 
necessary and forced upon us by the inner logic of Theory as we understand it in its 
essences -rather than a patchwork hypothetical "fix" of its symptoms- should guide 
our willingness to entertain new models: given that the "enveloping theory space " 
of established Theory is almost by definition multiply infinite. If we seek to avoid 
the (in our view very fortunately found) tight corners into which we are successful 
in driving the simplest -otherwise viable - model by changing the model itself so 
completely that the entire previous simplifying and constraining understanding must 
be jettisoned then we foreclose any opportunity of hearing the "tiny voice within" (the 
theory that is !) by which nature may seek to inform us of the path to the decoding 
of these ultimate enigmas. 

The above pontifications may seem too sententious, pretentious or subjective, 
since one theorist's minimality is another theorist's monstrous pedantry! So we put 
it in more concrete terms : changing the MSGUT's AM Higgs structure by adding 
additional 54's or other Higgs multiplets will likely completely destroy the GUT scale 
simplicity and calculability achieved - after tedious labour - in the current formulation 
of the0 EH IHl El E| MSGUT. We think it is preferable to first allow the 120-plet 
arbitrarily excluded by the BM-scenario (justifiably so at that early stage of searching 
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for the simplest viable model but now not an assumption worth protecting by ad hoc 
measures ): which leaves the GUT scale SSB unchanged. In this way one may hope 
to arrive at definite statements and falsifiability about the simplest viable model so 
as to define the next necessary (and falsifiable !) one. On the other hand extending 
the GUT scale structure to avoid the conflict with the FM data merely extends the 
catalog of non-falsifiable models without adding insight into the inner logic of the 
minimal theoretical structures adequate to encompass the data. 

Indeed, as shown in U the constrained analysis in the framework of the MSGUT 
points to the root of the problem lying in the multiple contradictory constraints 
faced by the 126 Yukawa couplings. A very simple calculable alternative, namely 
10 — 120 domination of the charged fermion sector coupled with very weak 126 
couplings that accentuate the Type I seesaw- due to its inverse dependence on these 
couplings through the large mass of the integrated out heavy right handed neutrinos - 
emerges as an obvious and well motivated resolution of the problem. Surprisingly just 
this possibility seems not to have been used to simplify the prima facie intractable 
nature of the fitting problem involving both the two symmetric Yukawa matrices of 
the 10, 120 and the antisymmetric couplings of the 120-plet. Indeed, guided by this 
insight, analysis of the scenario quickly yields rather novel insights into the structure of 
the fermion mass hierarchy. In particular, as shown by us in this paper,the mysterious 
maximality {Sin 2 9^ ~ 1.02 ± .04 ^7j) of the atmospheric mixing angle- which has 
no naturalness even in the BM Type I and Type II seesaw fits^ El El Ej ~~ here 
emerges as a direct consequence of the most accurately known fermion mass data at 
the GUT scale and the well motivated structural pattern chosen for the theory due 
to its previous debacle pQ ! 

Usually such "coincidences " are explained by means of somewhat thinly motivated 
additional discrete symmetries or textures. Yet here Nature seems once again to 
hint that the solution is both simple and "commonsensical" and at the same time 
more deep and profound than anything arbitrary speculation might have rigged. The 
fermion spectra at Mx seem seem to be fully compatible with the observed maximality 
of mixing without any assumption besides allowing all SO(10) FM Higgs types and a 
pattern of dominance (and "division of labour" ) that assigns each Higgs a consistently 
bearable "workload" : all serendipitously arranged so as to robustly yield maximal 
atmospheric mixing angle fixed at a nearly "geometric " value! 

Moreover the analysis has proved to be fully compatible with and robust under a 
systematic perturbation theory in which the CKM angle 8^ is considered as a common 
small number for structuring magnitudes of masses and angles (i.e universal Wolfen- 
stein parameter) in the whole hierarchy, leaving behind magnitudes 0(1) or smaller to 
be determined order by order in the expansion of the fitting problem. An expansion 
of of the three generation CP conserving case in this small parameter yields [21] a so- 
lution of the fitting problem that is based upon and respects the 2-generation analytic 
solution found in this paper as its robust core. Besides the remarkable b — r = s — \i 
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unification constraint-which strengthens the rough b — r equality long canonical in 
SUSY GUTs into the viable yet imminently falsifiable (i.e pending refinements in the 
bottom quark mass measurement etc) 6 — r pinning described above we find at 0(e 3 ) 
an additional constraint between angles : 6^ s = #12^23 > which is compatible with ob- 
servation at least as regards order of magnitude. Moreover since it arises at 0(e 3 ), 
by which order we expect CKM CP violation to rear its head, there is every reason 
to hope that this constraint will be modified to a more viable one involving the CKM 
phase when the full 3 generation analysis is completed |33j . Once the formalism of 
generic fitting problem has been clarified it will be ready for application in the Nu 
(realistic) MSGUT (NMSGUT) whose mass spectrum, zero mode Higgs couplings, 
and Baryon violating operators have been calculated [IB] within the framework and 
conventions of our decomposition [5] of SO(10) group theory to the maximal unitary 
sub group SU{4) x SU{2) L x SU{2) R (i.e the Pati-Salam group[3U|). 

Finally we note that the very nature of our ansatz dictates that the righthanded 
neutrino masses lie in the range below 10 12 GeV due to the smallness of the 126 cou- 
plings. Combined with the lower bound [25J that arises in the Leptogenisis models 
[2*o] naturally associated with SO(10) and the tight yet achievable pinning of rrib(Mx) 
to the more accurately known m T (Mx), we may hope that our approach has added 
two significant stable elements to the slowly but surely emerging picture of the in- 
tricate interconnections and tightly interdependent structures concealed withinin the 
Fermion mass Hierarchy. 

We earlier likened [Tj the circular linkages and mutual balance between the Large, 
Small and Geometric Mean mass scales to the fabled cosmic serpent that swallows 
its own tail : an ancient hermetic symbol expressing the same sentiment of wonder 
at the "boundary less" balance and intricate self sufficiency of the cosmic order that 
we seek to decipher. The wondrously patterned and enigmatic Fermion Mass pattern 
is then perhaps the supremely fascinating and valuable diadem that the ov p(36v poq 
1 is fabled to wear. What was still only the 'spoor of a grail' [H] now gleams -to our 
hopeful eyes - with the iridescence of the ouroborotic nagachudamani 2 which must -of 
course - also be our true philosopher's stone !. 
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8 Note Added 



On receiving an advance version of this paper, B. Bajc informed me that in their 
latest work on the 10 — 120 system including complex couplings for the 2 generation 
case they have also found restrictions on m& — m T very similar to those that form one 
of the central results of this paper [34]. 
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